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The hierarchy of moment equations derived from the nonlinear Boltzmann equa- 
tion is solved for a gas of Maxwell molecules undergoing a stationary Poiseuille flow 
induced by an external force in a pipe. The solution is obtained as a perturba- 
tion expansion in powers of the force (through third order). A critical comparison 
is done between the Navier-Stokes theory and the predictions obtained from the 
Boltzmann equation for the profiles of the hydrodynamic quantities and their fluxes. 
The Navier-Stokes description fails to first order and, especially, to second order in 
the force. Thus, the hydrostatic pressure is not uniform, the temperature profile ex- 
hibits a non-monotonic behavior, a longitudinal component of the flux exists in the 
absence of longitudinal thermal gradient, and normal stress differences are present. 
On the other hand, comparison with the Bhatnagar-Gross-Krook model kinetic 
equation shows that the latter is able to capture the correct functional dependence 
of the fields, although the numerical values of the coefficients are in general between 
0.38 and 1.38 times the Boltzmann values. A short comparison with the results 
corresponding to the planar Poiseuille flow is also carried out. 
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1 Introduction 



The Poiseuille flow, first studied by Hagen and Poiseuille towards the half of 
the 19th century, is still one of the classical examples in fluid dynamics [1]. It 
describes the steady flow in a slab or in a pipe under the action of a longi- 
tudinal pressure gradient, what produces a longitudinal macroscopic velocity 
with a typical parabolic profile in the transverse directions. Essentially the 
same effect is generated when the longitudinal pressure difference is replaced 
by a uniform external force mg (e.g. gravity) directed longitudinally. This lat- 
ter mechanism for driving the Poiseuille flow does not produce longitudinal 
gradients and so is more convenient than the former in computer simulations 
as well as from the theoretical point of view, especially to assess the reliabil- 
ity of the continuum description [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]. The first 
study of the Poiseuille flow driven by an external force we are aware of was 
carried out by Kadanoff et al. [2] , who simulated the flow with the FHP lattice 
gas automaton [17] to confirm the validity of a hydrodynamic description for 
lattice gas automata. In the context of a dilute gas, Esposito et al. [3] studied 
the Boltzmann equation and found that if the force is sufficiently weak there 
is a solution which converges, in the hydrodynamic limit, to the local equi- 
librium distribution with parameters given by the stationary solution of the 
Navier-Stokes (NS) equations. A generalized Navier-Stokes theory was seen 
to give a reasonable account of a fluid composed of molecules with spin when 
compared with molecular dynamics simulations [4]. 

The first kinetic theory analysis of the Poiseuille flow clearly exhibiting non- 
Newtonian behavior was carried out by Alaoui and Santos [5], who found an 
exact solution of the Bhatnagar-Gross-Krook (BGK) model kinetic equation 
for a particular value of g. The general solution of the BGK model under the 
form of an expansion in powers of g through fifth order was obtained by Tij and 
Santos [6]. The most interesting outcome of the solution, as first recognized by 
Malek Mansour et al. [7] , was the presence of a positive quadratic term in the 
temperature profile to second order in g, in addition to the negative quartic 
term predicted by the NS description. As a consequence of this new term, the 
temperature profile exhibits a bimodal shape, namely a local minimum at the 
middle of the channel surrounded by two symmetric maxima at a distance of 
a few mean free paths. In contrast, the continuum hydrodynamic equations 
predict a temperature profile with a (flat) maximum at the middle. The Fourier 
law is dramatically violated since in the slab enclosed by the two maxima the 
transverse component of the heat flux is parallel (rather than anti-parallel) to 
the thermal gradient. This phenomenon is not an artifact of the BGK model 
since the same results, except for the numerical values of the coefficients, 
were derived from an exact solution of the Boltzmann equation for Maxwell 
molecules through second order [8], as well as from approximate solutions of 
the Boltzmann equation for hard spheres by Grad's moment method [9,10]. It 
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is interesting to note that this correction to the NS temperature profile is not 
captured by the next hydro dynamic description, namely the Burnett equations 
[11]. Therefore, in order to describe a non-monotonic temperature field [which 
is an 0(g 2 )-order effect] from a hydrodynamic description one needs to go 
at least to the super-Burnett equations [8]. A systematic asymptotic analysis 
of the BGK equation with diffuse boundary conditions to second order in 
the Knudsen number [12] validates the results of Ref. [6] as corresponding 
to the Taylor series expansion of the normal (or Hilbert) solution around the 
center of the gap. The theoretical predictions of a bimodal temperature profile 
have been confirmed at a qualitative and semi- quantitative level by numerical 
Monte Carlo simulations of the Boltzmann equation [7,11] and by molecular 
dynamics simulations [9,13]. In the case of a dense gas, even though the non- 
monotonic behavior of the temperature profile may disappear, the existence 
of a quadratic term coexisting with the classical quartic term is supported by 
molecular dynamics simulations [14]. It is worth mentioning that when the 
Poiseuille flow is driven by a longitudinal pressure gradient rather than by an 
external force, the NS description is in much better agreement with Monte 
Carlo simulations of the Boltzmann equation [15]. 

Practically all the studies about the Poiseuille flow driven by an external force 
have considered the planar geometry [2,3,4,5,6,7,8,9,10,11,12,13,14,15], i.e. the 
fluid is assumed to be enclosed between two parallel, infinite plates orthogonal 
to the y-direction, the force acting along the z-direction. In the stationary 
laminar flow the physical quantities have a unidirectional dependence on the 
y coordinate only. While the planar geometry is the simplest one to study the 
Poiseuille flow, it is much more realistic to consider that the fluid is enclosed 
in a cylindrical pipe, the external force being directed along the symmetry 
axis z. In that case, the quantities depend on the distance r = ^/x 1 + y 2 from 
the axis and the fluxes have in general both radial and tangential components. 
Two of us have recently solved the BGK model for the cylindrical Poiseuille 
flow through fourth order in g [16]. The results showed that the structure of 
the hydrodynamic and flux profiles in the pipe geometry is similar to that of 
the slab geometry. In particular, the temperature exhibits a non-monotonic 
behavior as one moves apart from the pipe axis. On the other hand, the results 
are quantitatively sensitive to the geometry of the problem. For instance, 
the normal stress along the gradient direction is uniform in the planar case 
(P yy = const), while it has a non-trivial spatial dependence in the cylindrical 
case (P rr 7^ const). Also, comparison with the BGK solution for the planar 
geometry [6] shows that the differences between the kinetic theory results and 
the NS predictions are in general more pronounced in the planar case than in 
the cylindrical case. 

As is well known, in the BGK model kinetic equation the complicated non- 
linear structure of the Boltzmann collision operator is replaced by a single 
relaxation-time term to the local equilibrium distribution [18]. This is usually 
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sufficient to account for many of the nonequilibrium properties of the under- 
lying Boltzmann equation, at least at a qualitative level [19]. However, the ex- 
istence of only one effective collision frequency does not allow the BGK model 
to describe quantitatively those states where momentum and heat fluxes co- 
exist and are inextricably intertwined, as happens in the Poiseuille flow. This 
limitation of the BGK model is responsible for an incorrect value Pr = 1 of the 
Prandtl number [18], in contrast to the Boltzmann value [18,20] Pr ~ |. In the 
planar Poiseuille flow, it is possible to assess the validity of the BGK solution 
[6] by comparing it with the results derived from the Boltzmann equation for 
Maxwell molecules [8] and hard spheres [9,10]. The aim of this paper is to 
fill the gap caused by the absence of results from the Boltzmann equation for 
the cylindrical Poiseuille flow. We consider a gas of Maxwell molecules and 
solve the infinite hierarchy of moment equations stemming from the Boltz- 
mann equation through third order in the external force. The results confirm 
the correctness of the functional dependence of the hydrodynamic and flux 
profiles obtained from the BGK model [16]. On the other hand, the numerical 
values of the coefficients are in general different in both kinetic theories, as 
expected. The organization of the paper is the following. Section 2 is devoted 
to the description of the flow under study and its solution in a hydrodynamic 
description to NS order. The kinetic theory description is presented in Sec. 3. 
In order to solve the moment hierarchy in a systematic and recursive way, we 
carry out a perturbation expansion in powers of the external field in Sec. 4, the 
technical details being relegated to Appendix B. The results are extensively 
discussed and compared with those of the NS and BGK descriptions in Sec. 
5. Finally, the main conclusions of the paper are briefly presented in Sec. 6. 



2 Navier— Stokes description of the cylindrical Poiseuille flow 

The cylindrical Poiseuille flow studied in this paper refers to a monatomic 
dilute gas enclosed in a long pipe of radius R. Let us take the z-axis as the 
symmetry axis of the cylinder. The gas particles are subjected to the action 
of a constant external force per unit mass g = gz (e.g. gravity) parallel to 
z. After a certain transient stage, the system reaches a steady laminar flow. 
This nonequilibrium state is characterized by gradients of the hydrodynamic 
variables along the directions x and y orthogonal to the cylinder axis. By 
symmetry, the hydrodynamic fields are expected to depend on x and y through 
the distance r = ^/x 2 + y 2 from the z-axis. In our analysis, we are interested 
in the bulk region of the flow. This means that the radius R of the cylinder is 
assumed to be large enough (as compared with the mean free path) to allow 
for the existence of such a region. 

The steady- state balance equations expressing the conservation of momentum 
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and energy are 



V • P = pg, 



(1) 



V-q+P:Vu = 0, (2) 

where p is the mass density, u = u z z is the flow velocity P is the pressure 
tensor and q is the heat flux vector. The geometry of the problem suggests 
the use of cylindrical coordinates (see Appendix A). Using Eqs. (A.8)-(A.12), 
Eqs. (1) and (2) yield 

|- (rP rr ) = P H , (3) 
l£(rP r ,)=pg, (4) 

P ^ + I M = °' (5) 
Equations (3)-(5) are exact, but they do not constitute a closed set. 

In the Navier-Stokes (NS) description, the momentum and heat fluxes are 
assumed to be linear functions of the hydrodynamic gradients [1,18,20]. In the 
geometry of the cylindrical Poiseuille flow the NS constitutive equations are 
written as 

Prr = P<j><j> = P zz = P, (6) 

TZ = ' (7) 

dT 

Qr = -«^T, (8) 

q z = 0. (9) 
In Eq. (6), p — |trP is the hydrostatic pressure. Equation (7) is Newton's 
friction law, p being the shear viscosity. In addition, Eq. (8) is Fourier's law, 
where k is the thermal conductivity and T is the temperature. The latter is 
related to the number density n = p/m (where m is the mass of a particle) 
and the pressure p by the (local) equilibrium equation of state. In particular, 
for a dilute gas, p = nksT, ks being the Boltzmann constant. Equation (6) 
implies the absence of normal stress differences in a sheared Newtonian fluid, 
while Eq. (9) means that the heat flux is parallel to the thermal gradient 
in a fluid obeying Fourier's law. The combination of Eqs. (3)-(8) gives the 
following closed set of hydrodynamic equations: 
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In principle, the NS hydrodynamic equations (10)-(12) are too complicated 
to obtain its explicit solution for arbitrary g. On the other hand, the solution 
can be found as a series expansion in powers of g. To third order the result is 

p(r)=p , (13) 



^>=*°-^ i+ iysr i W ) ' (i4 » 




2 r ,2 



-±^r* + 0(g% (15) 
64 ?7o k 

where we have particularized to a dilute gas of Maxwell molecules, in which 
case rj oc T and k oc T [18,20]. The subscript in Eqs. (13)— (15) denotes 

quantities evaluated at r = 0, i.e. u = u z (0), T = T(0), r/o = rj(T ) : The 

corresponding expressions for the fluxes are 

+ (16, 

*M = i^r 3 + 0(9*), (17) 
16 r] 

complemented with Eqs. (6) and (9). It is interesting to note that the spatial 
variable r can be eliminated between Eqs. (14) and (15) to obtain the following 
nonequilibrium "equation of state" : 

T = T -^(u z -u ) 2 + O(g A ). (18) 
4k 

By equation of state we mean in this context a relationship holding locally 
among the hydrodynamic fields (p, u z and T) and that does not contain g 
explicitly (at least up to a certain order). Since the pressure is uniform in the 
NS description, it does not enter into Eq. (18). Equation (18) shows that in 
the Poiseuille flow the velocity and temperature profiles are closely related. 

From Eqs. (13)— (15) we can obtain the mass rate of flow in the NS description. 
Let us consider a circular section of radius a centered on the cylinder axis and 
orthogonal to it. The mass of fluid flowing across this surface per unit time is 
defined, in the reference frame moving with the fluid at r = a, as 

M(a)=2vr /° dr rpir) \u z (r) - uJa)] . (19) 
Jo 
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Using Eqs. (13)-(15) and taking into account that p = mp/k B T, one gets 

Before closing this Section, let us make a few remarks about our use of the 
term "Navier-Stokes solution" . By that we mean the solution to the coupled 
set of hydrodynamic equations (10)-(12). In fluid mechanics, on the other 
hand, it is usual to adopt a more restrictive point of view in which by Navier- 
Stokes equation one means the momentum equation with a constant viscosity 
77 and a constant density p ("incompressible fluid"). In that case, the solution 
of Eq. (11) would simply be 

u z (r) =u - ^r 2 , (21) 

regardless of the value of g. Accordingly, the mass rate of flow would be 

M(a) = (22) 

Equation (21) gives the typical parabolic velocity profile under Poiseuille flow, 
while Eq. (22) is known as Poiseuille's law [1]. Nevertheless, the shear viscos- 
ity i] depends in general on the temperature T and so the assumption of a 
constant viscosity implicitly implies a constant temperature. The same con- 
clusion follows from Eq. (13) and the incompressibility assumption p = const. 
The energy equation (12), however, shows that a velocity profile is compatible 
with a constant temperature only if the thermal conductivity k is allowed to 
be infinite. As a matter of fact, Eqs. (14) and (20) reduce to Eqs. (21) and 
(22), respectively, if we formally set k — > 00. Since in a dilute gas the relative 
strength of the thermal conductivity is comparable to that of the shear vis- 
cosity (as measured by the Prandtl number Pr = 5/cs?7/2m/t ~ |), the energy 
equation (12) cannot be neglected and hence T 7^ const. 

In summary, in the remainder of this paper we will refer to Eqs. (10)-(12) 
as Navier-Stokes equations and will compare their solution for a dilute gas 
of Maxwell molecules [Eqs. (13)-(15)], along with the momentum and energy 
fluxes [Eqs. (6), 9), (16) and (17)], with the results derived from the Boltzmann 
equation. 



3 Kinetic theory description 

In kinetic theory, the relevant information about the nonequilibrium state of a 
dilute gas is contained in the one-body distribution function /(r, v; t) [20]. Its 
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temporal evolution is governed by the non-linear Boltzmann equation [18,20], 
which in standard notation reads 



^/ + v-V/ + g-— f = fdv 1 f dn|v-vi|a(|v-vi|,fl)[/7f-//i] 

= J[f,f]- (23) 

The influence of the interaction potential appears through the dependence 
of the differential cross section a on the relative velocity |v — Vi| and the 
scattering angle 9. In particular, for Maxwell molecules cr(|v — Vi|,0) oc |v — 
vil" 1 . 

The knowledge of / in a dilute gas is sufficient to determine its macroscopic 
state. For example, the hydrodynamic variables and their fluxes are just ve- 
locity moments of the distribution function. Thus, the local number density n 
is defined by 

n = Jdwf. (24) 

Analogously, the flow velocity u and the hydrostatic pressure p are given by 
the following expressions: 

u = -fdwf, (25) 
nJ 

p=^JdvV 2 f, (26) 

where the peculiar velocity V = v — u has been introduced as the velocity 
of a particle relative to the flow velocity. The momentum and energy fluxes 
are given by the pressure tensor P and the heat flux vector q, respectively, as 
follows: 

p = mJdvVVf, (27) 



m 
~2 



JdvV 2 Vf. (28) 



In the case of the stationary Poiseuille flow in a pipe, the Boltzmann equation 
(23) becomes 

V ^ +V ^ +9 d^ z =J[fJ] - ^ 
This equation is invariant under the transformations 

(x, v x ) < — >■ (-x, -v x ), (y, v y ) < — >• (-y, -v y ), (v g , g) < — »• (-v z , -g). 

(30) 

Strictly speaking, Eq. (29) must be supplemented by the appropriate bound- 
ary conditions at r = R describing the interaction of the particles with the 
cylinder surface. However, in this work we are interested in the bulk region 
of the system (0 < r < R — 5, where 5 is the width of the boundary layer), 
which is not affected by the details of the boundary conditions. This implies 
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that the radius R of the cylinder must be sufficiently large, compared to the 
mean free path of a particle, to allow the existence of such a region (small 
Knudsen number). Here we will assume that this is the case and will look 
for the Hilbert-class or normal solution to Eq. (29), namely a solution where 
all the spatial dependence of the distribution function takes place through a 
functional dependence of / on the hydrodynamic fields n, u and T. 

Since the hydrodynamic variables of the gas and the associated fluxes are the 
first few moments of the function /, it is convenient to consider the hierarchy 
of moment equations stemming from the Boltzmann equation (29). A moment 
of an arbitrary degree k = k\ + k 2 + k 3 is defined as 

M klMM (x,y) = JdvV x k X 2 V z k3 f(x,y : v). (31) 

The invariance properties (30) imply the parity relations 



M klMM (x,y;g) = (-l) 1 M klMM (-x, y; g) 
= (-l) k2 M klMM (x, -y-g) 

= (-l) k *M klM , k3 (x,y--g). (32) 

Although the problem at hand calls for the use of cylindrical coordinates, we 
will use Cartesian coordinates to construct the moment hierarchy. The reason 
is two-fold. First, the collision integrals are easier to express in Cartesian (or 
even spherical) coordinates than in cylindrical coordinates. Second, the use of 
Cartesian coordinates will serve us as a test of the calculations since, when 
obtaining the cylindrical components according to Eqs. (A.l) and (A. 2), the 
relevant quantities must depend on x and y through the variable r = ^/x 2 + y 2 . 

The Boltzmann hierarchy of moments is obtained by multiplying both sides 
of Eq. (29) by V kl Vj 2 V ks and then integrating over the velocity space. The 
resulting expression is 



d d du z du z 

—M kl+ i >k2)k:i + —M kuk2+ i )k:i + fc 3 -^— M fcl+ i ife2ifc3 _i + fc 3 -^— M fcl ,fc 2 +i,fe3-i 

- gk 3 M kl;k2tk3 ^i = Jfci,fc 2 ,fc 3 ; (33) 

where 

Jkmm = J dW x k ^V z ^J[f, /]. (34) 

In the sequel, we will use the roman boldface k to denote the triad {k±, k 2 , k 3 } 
and the italic lightface k to denote the sum k\ + k 2 + k 3 . Thus, M k = M kl ^ 2tk3 
is a moment of degree k = k\ + k 2 + k 3 . 

The collision integral Jk involves in general all the velocity moments of the 
distribution (including those of a degree higher than k) and its explicit ex- 
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pression in terms of those moments is unknown. An important exception is 
provided by the Maxwell interaction potential tp(r) = K/r 4 . In that case, the 
collision rate is independent of the velocity, i.e. |v — vi|er(|v — Vi|, 6) = <J O (0), 
and thus J k can be expressed as a bilinear combination of moments of degree 
equal to or smaller than k [21,22]: 



J* = E C* k „M k ,M k „, (35) 



t^-fk 
k',k" 

where the dagger denotes the constraint k! + k" = k. The coefficients C k / k /; 
are linear combinations of the eigenvalues [22,23] 



\ r £ = J dVt a (6) 



1 + 5 r0 5 m - cos 2 ^ 9 -P e (cos 6 -j - sin 2 ^ 9 -P t (sin °- 

'(36) 

of the linearized collision operator, where Pg(x) are Legendre polynomials. The 
thermal conductivity and shear viscosity for Maxwell molecules (first obtained 
by Maxwell himself) are [20] 

5k B v 

K \T) = ; — , (37) 

v 1 2mn\ n 

V(T) = -f-, (38) 

where A 2 = |An = 0.436 x 3n^j2K/m. The eigenvalues \ r £ can also be 
expressed as linear combinations of the integrals [21] 



[ dna (e)sm 2s ^cos 2s ^-. (39) 



i2s 



In particular, A 02 = 3A 2 , A 04 = 7A 2 - fA 4 and A 06 = ™A 2 - ^-A A + ^Aq. 
The first few ratios A 2s /A 2 are A 4 /A 2 = 0.15778 . . ., A 6 /A 2 = 0.031196 ... and 
As/A 2 = 0.0066540. . .. The general expression (35) for J k was elaborated by 
Truesdell and Muncaster [21], who also give J k up to k = 5. A simplification 
algorithm for this formula, along with explicit expressions of J k for k = 6-8, 
has been recently proposed [24]. 

To simplify the subsequent analysis, it is suitable to use dimensionless quan- 
tities. First, let us define an effective collision frequency as 

v = n\ 02 . (40) 

Next, without loss of generality, we choose a reference frame stationary with 
the flow at r = 0, i.e. u = 0. Now we introduce the dimensionless variables 



* — 1 — 1 * — 1 * — 1 T~>* rji — 1 m 

9 = v u Q g, u z = v u z , p =p p, T = T T, 

M* = n -\- fe M k , J*=n W<4, r* = u v Q l v. (41) 
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In the above expressions, Vq = (kBTo/m) 1 / 2 is the thermal velocity and, as 
in Eqs. (13)-(18), the subscript means that the quantities are evaluated at 
r = 0. The reduced (gravity) acceleration g* has a clear physical meaning. 
The quantity h = v\j g is the so-called scale height [25], i.e. the characteristic 
distance associated with the external (gravity) field. Thus, g* = (v /vo)/ho 
represents the ratio between the mean free path and the distance through 
which a typical particle undergoes the action of g. While the parameter g* is 
a measure of the field strength on the scale of the mean free path, the Froude 
number 

Fr=^_ = ^ (42) 
gR R K J 

measures the strength on the scale of the radius of the pipe. The relationship 

between both dimensionless parameters is 

FY = «±, (43) 
9* 

where Kn = (vq/uq)/R is the Knudsen number. 

The solution of the hierarchy of equations (33) is a very difficult task, even for 
Maxwell molecules, due to its nonlinear character and the fact that moments of 
degree h are coupled to the first spatial derivatives of moments of degree k + 1. 
On the other hand, we will restrict ourselves to cases where the parameter g* 
is weak. This assumption is well justified in most of the practical situations, 
where values g* > 10~ 2 would be abnormally large. For instance, in the case 
of air at room temperature under the action of terrestrial gravity, g* ~ 10~ 12 . 
Consequently, we can solve the hierarchy (33) by means of a perturbation 
expansion in powers of g*, in a similar way as done in other works [6,8,16,26]. 



4 Expansion in powers of the external force 



For the sake of clarity, the asterisks will be dropped in this Section, so all the 
quantities will be understood to be expressed in reduced units [cf. Eqs. (41)], 
unless stated otherwise. 

Before performing the power series expansion of the moments, it is convenient 
to take into account the symmetry properties (32). For example, it is obvious 
that the profiles of the temperature and the hydrostatic pressure are even 
functions of g, while the flow velocity u z is an odd function. Therefore, we can 
write 

p = +P ( V + -" , (44) 

T = 1+T(V + ^ ( V + --- , (45) 
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u z = u m g + u i3) g 3 + ■■■ . (46) 

The expansion of the number density n = Mqqq = p/T can be obtained from 
Eqs. (44) and (45). These hydrodynamic fields are (even) functions of x and 
y through the distance r = ^Jx 2 + y 2 - As for the moments Mk, they are even 
(odd) functions of g if their index k 3 is even (odd). Thus, 

M k = M£ 0) + Mi 1] g + M^g 2 + g 3 + ■■■ , (47) 

where = if s + k 3 = odd. The zeroth-order term is the moment 

of the equilibrium distribution function normalized to p^ = 1, = 1. Its 
expression is 

Mf = C kl C k2 C k3 , (48) 

where C k = (k - 1)!! = (k - 1) x (k - 3) x (k - 5) x • • • x 1 if k — even, being 
if k = odd. 

In the first stage of the calculations, the expansions (44)-(47) are inserted into 
Eq. (33). By equating the terms of degree s in g on both sides, we get 



K*-l)/2] 

E 

i=0 



dx fc i+ 1 . fc 2,fc 3 -i 



+ ^ M (,-i-», 



(49) 



where [a] means the integer part of a. According to Eq. (35), the right side 
of (49) is 



k'.k" i=0 



t^-yk 



(50) 



Equation (50) is still very difficult to solve in general. First, it is nonlinear 
and couples moments of degree k and lower to moments of degree k + 1. 
Second, the coefficients are unknown functions of the spatial variables 
x and y. To overcome the latter obstacle, we assume that a self-consistent 
solution of Eq. (49) exists in which the coefficients in the series defined by 
(44)-(47) are polynomials in x and y. This ansatz is justified by the results 
obtained in the case of the planar geometry from the BGK model [6] and the 
Boltzmann equation [8], as well as in the case of the cylindrical geometry from 
the BGK model [16]. More specifically, we assume that the coefficients of the 
hydrodynamic profiles are of the form 

P (s) =I>2?r 2 \ (51) 
i=i 
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T^ = ±T^r 2 \ (52) 
i=i 

« W = E4f^. (53) 
i=i 

The task becomes more difficult when guessing the spatial dependence of the 
coefficients It can be checked that the only representation that leads to 

consistent solutions of the hierarchy (49) is of the form 

j=0 j=0 

where a(s = 1) = 1 and a(s) = 2s for s > 1. From Eqs. (32) it follows that 
XkiMM = if s + fc 3 = odd or i + fci = odd or j + k 2 = odd. Insertion of Eq. 
(54) into Eq. (50) shows that has the same structure as namely 

4 S) = E E -C "V V- (55) 

i=0 j=0 

The numerical coefficients ^2i\ u 2i an d Xk' 1 '^ are determined by con- 
sistency. To do so, we insert Eqs. (51)-(54) into Eq. (49) and equate the 
coefficients of the terms x l yi in both sides. This yields a hierarchy of linear 
equations for the unknowns. In the hierarchy, the coefficients of the moments 
of degree k + 1 , k and k — 1 on the left-hand side are related to those of degree 
k, k — 1, k — 2, ... on the right-hand side. In addition, the coefficients of the 
type X^+~ik2k 3 an d xii'fc2+i^3 011 ^ ne left-hand side determine those of the 
type XkiMM on ^ e r ight-hand side. In spite of the intricacy of the resulting 
hierarchy, it can be recursively solved by following steps similar to those made 
in the case of the planar geometry [8] . To get explicit results to order s = 3 it 
is necessary to consider the moment equations up to degree k — 8. 

The method is outlined in Appendix B and the main results are presented in 
the next Section. 



5 Discussion 

5.1 Hydrodynamic profiles 

By proceeding along the recursive scheme described in Appendix B, we have 
obtained the explicit spatial dependence of the coefficients in the series ex- 
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pansions (44)-(46) through order g A and those of the series expansion (47) for 
moments of second (k = 2) and third (k = 3) degrees through order g 3 . Here 
we display the results to third order. Also, to ease the physical interpretation 
of the results, they are expressed in cylindrical coordinates and in real units. 
The hydro dynamic profiles are 



Por 2 g 



4r/ 



l + ^^4r 4 + cJ^% 2 + C PO?7oV 



144 k t]qT 



kBT 



" ; pI 



+ 0(g 5 ), 
(56) 



p(r) = po 



1 + " 2 



+ 0(g% 



T{r) = T 



1 Po9 2 I . / "'.'/ . .2 



64 K r] T 



r 4 + Ci 



k B T 



+ 0(g% 



(57) 



(58) 



where ( u = ^§§ ~ 0.35, C = 24.322 . . , C P = ^ and Ct = ^ ^ 0.19. We 
recall that the subscript in a quantity represents its value on the z-axis, i.e. 
at r = 0. Elimination of r among Eqs. (56)-(58) allows one to get the following 
nonequilibrium equation of state: 



Vo 



Ct T 



T = T --^(u z - u Y + ^^(p- Po ) + Q(g A ). 



CpPo 



(59) 



The non-zero elements of the pressure tensor are 



P rr {r) = Po 



1 



1 + Q Cp \k B T 



-"i \ 2 2 /-i pmW 

r — SB 



Po 



+ 0{g% 



(60) 



P<M>(r) = Po 

Pzzij-) = Po 



7 fmg_ 



P 3 

Po . 

PovW 



r 2 + 2(> P , 
Po 



+ 0(g% 
+ 0(g% 



(61) 
(62) 



P n (r) = f r 



1 + 



1 PoV „4 , C P -Ct 



-r + 



mg 



192 k ?7oTo ' 2 V k B T 

where (' P = 1.7388 Finally, the components of the heat flux are 



+ 0(g% (63) 



2^2 



q r (r) 



1 Po9 

16 T]o 



■r 3 + 0(g 4 ), 



(64) 



Qz(r) = ~ 



2 Komgf 
5 kn 



Co-C 



Pofl 2 



r 2 -C 



, r p oVo9 2 
Po 



+ 0{g% 



(65) 
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where Co = 1, Cg = § - 0.29, (' q = 3.2527 ... and Q = 86.415 . . .. 

Equations (56)-(65) are the main results of this paper. They give the exact 
profiles in the bulk region of the hydrodynamic fields and their fluxes for a 
dilute gas of Maxwell molecules through third order in the external force. 
Comparison of Eqs. (6), (9) and (13)-(18) with Eqs. (56)-(65) show that the 
NS constitutive equations predict ( u = (' u = ( p = Cr — Cp — Co — C? = 
C = = ^' These vanishing values of the coefficients summarize the main 
limitations of the NS description of the Poiseuille flow induced by an external 
force. Equation (65) shows that q z ^ 0, even to first order in g (Co 7^ 0), 
despite the absence of a thermal gradient along the longitudinal direction. 
This represents an obvious breakdown of Fourier's law and is in fact a Burnett- 
order effect [8]. To second order in g the hydrostatic pressure is not uniform 
but increases radially (C P > 0); analogously, there are normal stress differences, 
so that P rr < P^ < p < P zz . These are non-Newtonian effects. 

Also to second order in g, there exists a positive quadratic term (Cr > 0) 
in addition to the negative quartic term in the temperature profile. The for- 
mer term is responsible for the fact that the temperature has a local min- 
imum rather than a maximum at r = 0. The maximum temperature is lo- 
cated at a distance r max = v / 32Ct^o — 2.49£o (independent of g), where 
£ = (rjoKoTo) 1 / 2 / po = y / |Pr(t> /t / o) is the (local) mean free path at r = 0; 
the relative difference between the maximum temperature and the tempera- 
ture on the cylinder axis is (T max - T )/T = lQ&(£ /h ) 2 ~ O.QOA(£ /h ) 2 , 
where we recall that h = kBT /mg. This result represents a dramatic viola- 
tion of Fourier's law: while the heat flows radially outwards, the temperature 
increases from r = to r = r max . Therefore, within the region < r < r max , 
the heat flows from the colder to the hotter layers. This paradoxical effect, 
which is beyond the Burnett description [11], does not violate the conserva- 
tion of energy (5) since the radial increase of rq r (which tends to cool the 
gas) is exactly compensated for by the viscous heating term rP rz du z /dr. It is 
interesting to note that from Eqs. (58) and (64) one can obtain 

dT 1 

~ K ~dr =Qr ~ 9 r -- V2 ^ + ° ( ^ 4) ' (66) 

where V 2 X = r~ l d(rdX/dr)/ dr is the Laplacian in cylindrical coordinates. 
Equation (66) is an extension of Fourier's law (8) showing that the sign of the 
thermal gradient results from a competition between the radial component of 
the heat flux and its Laplacian [10]. The latter dominates for r < r max and so 
dT/dr > in that region. 

Third-order contributions appear in u z and P rz but they do not have a qual- 
itatively important influence. They are responsible for a spatial variation of 
the flow velocity u z and the shear stress P rz more pronounced than expected 
from the NS equations. The same happens for the mass rate of flow defined 
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by Eq. (19). Insertion of the hydro dynamic profiles (56)-(58) yields 



M(a) = ^ [l + A^L a 4 + jtg + C ~ Cr (rng_\ 2 fl2 ^ Pjrtf 
y ' 8 Vo [ 38A Vo k T 3 \k B Tj Sn pi 

+0(g 5 ). (67) 

Comparison with Eq. (20) shows that the mass rate of flow grows with the 
radius a more rapidly than predicted by the NS approximation, but otherwise 
the qualitative behavior is similar in both descriptions. 



5.2 Generalized constitutive equations 



Let us use Eqs. (56)-(65) to rewrite the fluxes under the form of generalized 
constitutive equations. We begin with the shear stress: 



Prz = -V 



X 



du z 
dr 

£ 2 



1 - 



c: 



'2 4 

' V T — 



15( T T 

2l 



Cu 

15 I 15Cr 



+ 8C + 6Ct - 2C P 



du 



ksT /m \ dr 



+ 0(g 5 ). 



(68) 



Here i] oc T is the local shear viscosity at r, so rj — i] oc T — T = 0(g 2 ). 
Analogously, £ = {r]KT) 1 ^ 2 /p is the local mean free path and £ — £ = 0(g 2 ). 
In Eq. (68) the NS term is of order g, while the super-Burnett terms (i.e. 
of third order in the hydrodynamic gradients) are of order g 3 . The normal 
stresses can be written as 



P 



1 + 



3 T 



V 2 T 



15C P 3 J p 



V 2 p 



+ 0(g 4 ), (69) 



P<H> =P 



2rr ( Cp 



1 + 2C V — V Z T - , 

T \15( p 



P 



1 - 



16C P l 2 
3 T 



V 2 T + 



T | -V 2 p 



+ 0(g% 



, 15 C P 3 / p 



(70) 



+ 0(^ 4 ). (71) 



Thus, in order to obtain the normal stresses to order g 2 we need to include at 
least the Burnett contributions in the Chapman-Enskog expansion. 



The radial component of the heat flux can be expressed as 
ft- = -«|: (T + AC T £ 2 V 2 T)+0(g 4 ). 



(72) 
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In this case, both the NS term and the super-Burnett term are of order g 2 . As 
for the longitudinal heat flux, one has 




(73) 



where 6 3 = -3, A = 3( = 3, 5 = \{l + 64^) = f and 

fi = + IQCu + 2Ct(7 + 64C,) + 4C; - 8C P ^ 64.141 .... (74) 



The Burnett term headed by 64 is of of first order in g, while the Burnett terms 
with 63 and 8 5 are of third order. In addition, there are super-super-Burnett 
terms that are also of order g 3 . 

Equations (68)-(73) can be written in other equivalent forms by taking into 
account that the hydrodynamic gradients are not independent in our problem. 
For instance, 



V 2 P = -4C P ^ 2 V 4 T + Otf) = yX v ^ 2 (VV) 2 + 0{g% (76) 



In particular, Eq. (76) implies that the Burnett term V 2 p and the super- 
super-Burnett terms V 4 T and (V 2 w 2 ) 2 are of order g 2 . Therefore, in order to 
obtain correctly the profiles (56)-(65) from the Chapman-Enskog method one 
would actually need to go beyond the apparent orders in the gradients of Eqs. 



The previous analysis shows that the Chapman-Enskog expansion and the 
expansion in powers of g in the Poiseuille flow are quite different, even though 
the gradients are induced by the external force. Both expansions require a 
weak field, namely h ^> £ . On the other hand, while the Chapman-Enskog 
expansion requires that £qTq 1 X7 2 T 1, i.e. r ^> £ , the expansion in powers 
of g remains valid for distances smaller than or of the order of the mean free 
path. Since the non-monotonic behavior of the temperature profile occurs for 
< r < r max ~ £ , this effect is neglected by the NS description. 




) 



2 



t(!>-^ t ) +o(94) - (75) 



k B T \ dr 



(68)-(73). 



17 



Table 1 

Numerical values of the coefficients appearing in Eqs. (56)-(65), as obtained from 
the Boltzmann equation for Maxwell molecules (BM) and for from the BGK model 
for Maxwell molecules. The ratio between the BM and BGK values are displayed in 
the third column. 



c 


BM 


BGK 


BM/BGK 


U 


4429 
12600 


89 
200 


0.79 


Cm 


24.322 


64 


0.38 


Cp 


3 

To 


3 

TO 


1 


Ct 


34 
175 


7 

50 


1.38 


Cp 


1.7388 


92 
25 


0.47 


Co 


1 


1 


1 




37 
128 


15 
64 


1.23 


C' q 


3.2527 


209 
50 


0.78 


c; 


86.415 


5432 
25 


0.39 



5.3 Comparison with the BGK solution 

Once we have discussed the limits of validity of the NS theory (which otherwise 
provides a powerful tool to study the macroscopic state of many nonequilib- 
rium flows) in our Poiseuille problem, it is worthwhile wondering whether 
the BGK model kinetic equation succeeds in capturing the main qualitative 
features of the Boltzmann solution. The solution of the BGK model for the 
cylindrical Poiseuille flow was obtained by Tij and Santos [16] to fourth order 
in g for a general class of power-law repulsive potentials. The results agree 
exactly with the structure of Eqs. (56)-(65). On the other hand, most of the 
coefficients Cs have different numerical values. The failure of the BGK model 
to reproduce the exact Boltzmann values for those coefficients is not surpris- 
ing, given the simplicity of the model, which assumes that the practical effect 
of collisions is to make the distribution function relax to local equilibrium with 
a single characteristic rate v. In fact, it is well known [18] that the BGK model 
yields a value (Pr = 1) for the Prandtl number Pr = 5kBi]/2mK different from 
the correct one (Pr = |). Table 1 compares the numerical values of the coeffi- 
cients C's given by the Boltzmann equation for Maxwell molecules (BM) with 
those obtained from the BGK model. Except for C P and Co, the coefficients are 
different in both theories, the ratio being bounded between 0.38 and 1.38. For 
instance, the location r max of the maximum temperature is in the BGK model 
15% smaller than in the Boltzmann equation and the temperature difference 
T max — T in the BGK model is almost half the Boltzmann value. A comparison 
of the coefficients appearing in the hydrodynamic fields to fourth order in g 
[not shown in Eqs. (57) and (58), but given in Appendix B] shows that the 
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Fig. 1. Profile of the flow velocity u z (r) for g = Q.Ivq/Iq, as obtained from the 
Navier-Stokes (NS) equations (dotted line), the Boltzmann equation for Maxwell 
molecules (BM) (solid line) and the BGK model (dashed line). 

discrepancies between the BGK and Boltzmann equations increase with the 
power of g. In particular, the series obtained from the BGK model seem to 
diverge more quickly than those obtained from the Boltzmann equation. 

In order to illustrate graphically the differences between the NS, BGK and 
BM descriptions, we plot in Figs. 1-6 the hydrodynamic and flux profiles for 
the case g = 0.1vq/£ [which corresponds to ho/£ = 10], when only terms 
through third order in g are retained. Of course, higher order terms are not 
expected to be negligible for that particular high value of the external field 
[27], especially as one departs from the discussed in Subsection 5.5. 

However, the retained terms are sufficient to show the qualitative differences 
in the predictions of the three approaches, so that the main aim of Figs. 1-6 is 
to highlight those differences. It must be born in mind that the BGK curves in 
Figs. 1-6 correspond to Eqs. (56)-(65) with the BGK values for the coefficients 
C's (cf. Table 1) but with the correct value Pr = | of the Prandtl number. As 
expected, the BGK model is in qualitative agreement with the BM results. 

5.4 Comparison with the planar Poiseuille flow 

Let us now add a few comments on the similarities between the cylindrical 
and the planar Poiseuille flows. The solution of the Boltzmann equation for 
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Fig. 4. Same as in Fig. 1, but for the shear stress P rz (r). The NS and BGK 
are slightly below and above, respectively, the BM curve. 
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Fig. 6. Same as in Fig. 1, but for the components of the heat flux. Note that the 
NS, BM and BGK curves coincide for the radial component q r (r). 

Maxwell molecules to order g 2 in the planar geometry [8] yields profiles anal- 
ogous to those of Eqs. (56)-(65). In particular, the temperature profile is 

+ 0(g% (77) 

where now = 1.0153. As a consequence, the temperature has a maxi- 
mum with (T max - T )/T = 3C|(VM 2 - 3.09(VM 2 at \y\ = y max = 
V6Ct^o — 2.47£ - Therefore, the maximum temperature occurs at a separa- 
tion y max from the middle plane practically equal to the radial distance r max 
of the cylindrical case, but the relative temperature change (T max — To) /To 
is about five times larger in the former case than in the latter. Interest- 
ingly, this effect is accurately captured by the BGK model [16], as explained 
by the fact that in the planar Poiseuille flow one has [6,8] C BM /C BGK = 1, 
(Bm^bgk = 1.0153/(19/25) ~ 1.34, Cp BM /Cp BGK = 6.2602/(306/25) ~ 0.51, 
C( BM /C( BBGK = 1) m close agreement with the corresponding ratios in the cylin- 
drical Poiseuille flow (see Table 1). 



T(y) = To 



1 - 



12 K f?oTo 



y 



+ Ct( 



mg 



5.5 Boundary conditions 



When solving the Boltzmann hierarchy of moments [cf. Eq. (33)] by a per- 
turbation expansion in powers of gravity we have taken the hydrodynamic 
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quantities at the axis (r = 0) as reference values. In so doing, we have avoided 
the need of imposing specific boundary conditions. The price to be paid is that 
it is difficult to determine how far from the axis the profiles (56)-(65) remain 
valid. Consider for instance the temperature profile (58). The first neglected 
term 0{g 4 ) is a polynomial in r of degree 8 [cf. Eq. (B.19)]. This term can ac- 
tually be neglected versus the retained terms only if g* 4 (r/£ ) 8 <C g* 2 (r/£ ) 4 , 
i.e. if 

g*(r/£ ) 2 « 1. (78) 

The same conclusion is obtained from the velocity profile (56) and the pressure 
profile (57). If r ~ £q, condition (78) is equivalent to g* <C 1. On the other 
hand, a much stronger condition g* <C Kn 2 is needed to extend the profiles 
(56)-(65) to distances r comparable with the pipe radius R. In order to find 
the hydrodynamic fields for r ~ R when condition (78) does not hold, one 
would need to solve the full Boltzmann equation with the appropriate (e.g. 
diffuse reflection) boundary conditions at r — R corresponding to a given wall 
temperature T w . This has been done by Aoki et al. [12] in the case of the BGK 
equation for the planar geometry. They derived the hydrodynamic equations 
and the matching conditions corresponding to the normal (or Hilbert) solution 
by an asymptotic expansion in powers of the Knudsen number Kn and for 
values of the Froude number of order Fr ~ 1/Kn. In view of Eq. (43), this 
implies that g* ~ Kn 2 . The successive hydrodynamic equations had to be 
solved numerically, but Aoki et al. verified that a Taylor series expansion 
around the center of the gap agreed with the perturbation solution found in 
Ref. [6]. 

An asymptotic method analogous to that worked out by Aoki et al. [12] would 
be much more difficult to carry out in the case of the Boltzmann equation 
for the cylindrical geometry. Nevertheless, it is natural to expect that the 
results derived in this paper would be recovered by performing a Taylor series 
expansion around r = of the solution for g* ~ Kn 2 . This in addition would 
provide the matching conditions for the profiles (56)-(58) in the region r/£ ~ 
j/- 1 / 2 ~ KiT 1 . 



6 Conclusion 

In this work, we have studied the (laminar) stationary Poiseuille flow in a 
cylindrical pipe induced by an external force mg. We have considered a di- 
lute gas of Maxwell molecules and have analyzed the hierarchy of moment 
equations associated with the nonlinear Boltzmann equation. A consistent so- 
lution has been found under the form of a perturbation expansion in powers 
of g through third order. In principle, the method can be pushed to higher or- 
ders but, not only the algebraic intricacy of the solution scheme grows rapidly, 
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but the results would not be of practical use given the suspected asymptotic 
character of the series. The discussion of the results has focused on the most 
relevant physical quantities, namely the hydrodynamic profiles (flow velocity, 
pressure and temperature) and the associated fluxes (pressure tensor and heat 
flux). The results obtained here show once again that the Navier-Stokes (NS) 
theory remains unable to envisage the correct hydrodynamic profiles, even at 
a qualitative level, for distances smaller than or of the order of the mean free 
path. The most important limitation of the NS description is that it predicts 
a monotonically decreasing temperature as one moves apart from the cylin- 
der axis. In contrast, our solution to the Boltzmann equation shows that the 
temperature has a local minimum (T = To) at the axis (r = 0) and reaches a 
maximum value T = T max ~ To [1 + 0.6(£o/^o) 2 ] a t a distance from the center 
r = r max ~ 2.5£q of the order of the mean free path £ . In addition, a longitudi- 
nal component of the heat flux exists (q z ^ 0) in the absence of gradients along 
the longitudinal direction and normal stress differences (P rr < P^ < p < P zz ) 
are present. On the other hand, the BGK model provides results qualitatively 
similar to those found in this paper [16]. Notwithstanding this, the profiles 
predicted by the BGK model require some corrections from a quantitative 
point of view, as expected. In particular, the BGK model underestimates the 
maximum value T max of the temperature, as well as its location r max . From 
that point of view, one might say that the BGK model deviates from the NS 
predictions less than the Boltzmann equation, at least for Maxwell molecules. 
The same conclusion arises from a comparison between the Boltzmann and 
BGK results for the planar Poiseuille flow [6,8]. Nevertheless, it is worth not- 
ing that the BGK model is accurate in describing the effect of the geometry 
of the Poiseuille flow on the profiles. For instance, it correctly predicts that in 
the planar case the maximum temperature is located at a separation y max from 
the middle plane practically equal to the radial distance r max of the cylindrical 
case, although the relative temperature change (T max — T )/T is about five 
times larger in the former case than in the latter. 

Although we have not considered interactions different from Maxwell's, we 
expect the results found here to remain essentially valid, except that the nu- 
merical values of the coefficients may change. This expectation is supported 
by the results obtained with the BGK kinetic model [16]. In that case one 
finds that the coefficients ( p , (t, C'p an d Co are universal, while the coefficients 
C«, C Cg, C' q and Q change by 5%, 0.5%, 3%, 0.7% and 0.8%, respectively, 
when going from Maxwell molecules to hard spheres. 

Let us conclude by asserting that the stationary Poiseuille flow driven by an 
external force, both in a slab and in a pipe, is a good and conceptually simple 
example in fluid dynamics showing the limitations of a purely continuum the- 
ory in contrast to a kinetic theory approach. Next, we want to emphasize the 
important role played by the BGK model kinetic equation to pave the way 
towards a more fundamental theory based on the Boltzmann equation. 
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A Cylindrical coordinates 



The relationships between the cylindrical and Cartesian components of a vec- 
tor A and a tensor B are 



(A.1) 
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(A.2) 



U = -y/r x/r (A.3) 
, 1 
is a unitary matrix and is its transpose. 

In general, the divergence of a vector A in cylindrical coordinates is [28] 



___ . 19,., 19, 9 , 
V • A = -— (rA r ) + --^ + — A z . 
r or r 90 az 

The cylindrical components of the divergence of a tensor B are [28] 
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(A.6) 
(A.7) 
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In the Poiseuille problem the cylindrical components of the flow velocity u, 
the heat flux q and the pressure tensor P depend on the radial variable r only. 
Consequently, 

V.q = ~(nz r ), (A.8) 
(V • P) r = H (rP rr ) - ±.P„, (A.9) 
(V-P) = O, (A.10) 

(V • P), = ~ (rP„) , (A.ll) 

P:Vu = P„^. (A.12) 
In Eq. (A. 10) we have taken into account that, by symmetry, P r< ^ = P^ = 0. 



B Solution of the hierarchy of moment equations 



In this Appendix we outline the steps needed to solve the hierarchy (49) with 
the profiles (51)-(54). To do so, we proceed sequentially from a given order s 
to the next order s + 1. It is also necessary to use the consistency relation 

P=\ (M 2)0)0 + M 0)2i0 + M 0i o, 2 ) (B.l) 
to any order in g. Let us start from the first-order equations. 



B.l First order in g (s = 1) 



Setting s — 1, Eq. (49) becomes 



9 M (1) + 9 M (1) + k 9u(1) M (0) 

Q x m ki+1MM + o y m kiM+i,k 3 + K 3 q x M k 1+ i,k 2 ,k 3 -i 

Next, inserting (53) and (54) into (B.2), we obtain 
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(1,1,0) , v (i,o,i) , M (0) J^ + ot M (0) 7/ (1) 7/ 

(B.3) 

where J^' 1 '^ is defined by Eq. (55). Equating the coefficients of the same 
degree in x and y, Eq. (B.3) decouples into the following set of equations: 

= 4:i°l 3 , (B.4) 

a^ViA-rf = J £K 3 , ( B - 5 ) 

(1,1,0) (1,0,1) _ , M (0) _ ,(1,0,0) / R fi x 

A.fci+l,fc 2 ,fe3 ^ A.fci,fc 2 +l,fc 3 ft, 3- iw fci,fe 2 ,fe 3 -l — J ki,k 2 ,k 3 - \ D - U ) 

Let us start with the moments of small degree (k = 2). For example, by taking 
{ki, k 2 , h} = {1, 0, 1} in (B.4) and {h, k 2 , k 3 } = {0, 1, 1} in (B.5,) we obtain 

= -xKi 0) , (B.7) 

2«?> = - X ^f. (B.8) 
Now we insert the above equations into Eq. (B.6) with {hi, k 2 , k 3 } = {0, 0, 1}. 
This allows us to obtain the coefficient of the flow velocity profile to first order 
in g: 

«?> = -\- (B.9) 

The process of solution continues in the same manner for the coefficients of 
moments of degree k > 2. We use the relations (B.4) and (B.5) in a first stage 
to determine the coefficients Xk' 1 '^ an d Xk' '^ w hh k — 4. In the second 
stage, the obtained results are inserted into (B.6) to evaluate the coefficients 
of degree k — 3. The same steps can continue infinitely by following the routes 
A; = 6 — > A: = 5; A; = 8 — > A: = 7; .... It is not necessary to evaluate all the 
coefficients, but some of them are essential to the solution of the problem to 
the next order in g. 

B.2 Second order in g (s = 2) 

As done in the first-order evaluation, let us take s = 2 in Eq. (49) and insert 
Eqs. (51)-(54). We then obtain an equation more complicated than that of 
the first order. Of course, the symmetry property Xk * — 0iii + ki — odd or 
j + k 2 = odd must be applied. By equating the coefficients of the same degree 
in x and y we get the following set of algebraic equations: 
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a) i + j = 4 

b) i+j = 3 



= 4 2,i ' j) , (B.10) 
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e) i + j = 

(2,1,0) , v (2,0,l) _L V (1.0,0) _ r(2,0,0) (U-\A\ 

Xfci + l,fe 2 ,fe3 Afei,fc 2 +l,fe3 fi '3Afe 1 ,fe 2 ,ifc3-l — J k • l " 1 ^ 

In the above equations the convention Xk'' J ' = if i < or j < is implicitly 
understood. 

According to the system (B.11)-(B.14), the coefficients of degree k + 1 with 
i + j = a + 1 determine those of degree k with i + j — a. For instance, Eq. 
(B.10) allows one to determine the coefficients with k = 4 and % + j — 4. 
The results obtained are inserted into (B.ll) to determine the coefficients of 
degree k = 3 and % + j — 3. Then we go to Eq. (B.12) in order to determine 
those of degree k = 2 and % +j — 2. In general, the process of solution is done 
according to the following scheme: 

{xf +J = 4) } - {x£^= 3) } - {xK^ 2) } - {xK^} - {x£f^>}, 

(B.15) 

where {xi } denotes the set of coefficients of the same degree, namely 
{xi 2 ' Q ' ) } = {Xk ^i + ^2 + ^3 = k,i + j = a}. The expression of J^' 1 '^ 
is a linear combination of the coefficients Xk' !J \ ^4 2 \ ^i 2 "* an d P^- Once a 
row of the chain (B.15) is solved, it is necessary to go to the following row 
k + 2. This process continues until all the coefficients of the required profiles 
are completely determined. The results for the pressure and temperature are 

p^ = -r 2 , T< 2 > = — r 2 - — r 4 . (B.16) 
y 10 ' 175 240 V ; 



_ j(2,i,j) 
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B.3 Third and fourth orders in g (s = 3,4) 



The process to third order in g is practically the same as the previous one. Since 
the pressure and the temperature are even functions of g, only the coefficients 
of u z appear at this level. The result is 

M (3)_ (3) 2 4429 4 1 6 (Bl7) 

U ~ U2T 50 400 r 2160 r ' {BA7} 

where the coefficient = —6.0806 ... is related to the ratio A 4 /A 2 . 

The calculation algorithm is essentially the same for the higher orders in g. 
However, the equations become more and more cumbersome and their solution 
requires a considerable computational effort. Here we give the final results for 
the pressure and the temperature to fourth order: 

p W=p(V+pfV + ^ 5 r 8 , (B.18) 



T (4) = T W 2 + T (4) r 4 _ 30931 r 6 _ ^_ r 8 (B . 19) 

2 4 13 608 000 691200 v ; 

where pj, 4) = -24.160 . . ., pj 4) = -0.062666 . . ., T 2 (4) = -59.511 ... and T 4 (4) = 
—0.52856 . . .. The high values of \p^\ and {T^l suggest that the series (44)- 
(47) are only asymptotic, in agreement with the situation in the case of the 
BGK model [8]. 

Along with the hydrodynamic fields, we have obtained all the moments of 
degrees k — 3 and k = 4 to third order in g. This allows us to get the pressure 
tensor and the heat flux. Their expressions in cylindrical coordinates are given 
in Sec. 5. 
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